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Abstract 

A general framework for solving the Boltzmann equation for a 2-dimensional 
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(N . 

(^ , electron gas (2DEG) in random magnetic fields is presented, when the random 
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fields are included in the driving force. The formalism is applied to some 
recent experiments, and a possible extension to composite fermions at v = 1/2 
is discussed. 
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I. INTRODUCTION 

Considerable interest both theoretically and experimentally has recently focussed on 
the response of a 2DEG in spatially inhomogeneous external fields. Examples include the 
so-called Weiss-oscillations when a periodic electric field modulation is imposed,El similar 
effects are predicted for periodic magnetic fields,@ and random magnetic fields have been 
studied recently by several authors.ETEl In this paper we show how to set up a calculation 
of the transport coefficients of a 2DEG in a random magnetic field. The magnetic fields 
in the experiments cited above are all fairly low, well below the onset of Shubnikov-de 
Haas oscillations, leading us to expect a semiclassical treatment to be valid. We explicitly 
disregard any weak localization effects, which are well understood and can be subtracted 
from the experimental results. 

The plan of the paper is as follows: In section |H] we describe the basic physical system 



and our model. In section III we discuss how to include the random magnetic field in the 



Boltzmann equation. Section [IV] and [V| are devoted to the calculation of the resistivity 



tensor. Finally, in section [VT] we give some applications of our theory to recent experiments, 
and discuss suggestions for future experiments. We also briefly consider the applicability of 
our theory to the Quantum Hall state at v — 1/2. The conclusion is in section |V1J| . 

II. PRELIMINARIES 

Let us consider a 2DEG in a constant (external) electric field Eq, and a spatially varying 
magnetic field applied perpendicular to the plane of the electron gas: B{r) = (0, 0, B(r)) (r 
is in the plane). There have been different realisations of the magnetic field: One, originally 
proposed by Rammer and ShelankovQ for studying weak localization effects in inhomogeneous 
magnetic fields, consists in a type II superconducting gate on top of the heterostructure con- 
taining the 2DEG. When above the lower critical field B c \ for the superconductor, an applied 
field will penetrate in flux tubes, each containing an integral number of (superconductivity) 



flux quanta $q C = h/2e. Because of disorder in the superconductor the flux tube distribu- 
tion will not be the Abrikosov lattice, but a more or less random configuration, giving rise to 
a randomly modulated field in the 2DEG. Another approach was taken by Smith et alB who 
deposited small lead grains on the surface of a high mobility GaAs/AlGaAs heterostructure. 
For the grain sizes considered Pb is a type I superconductor. Below B c there will be (partial) 
flux expulsion from the grains, creating an inhomogenous magnetic field in the 2DEG. Still 
another possibility is to deposit small grains of a permanent magnetic material on top of 
the 2DEG. The different kinds of inhomogeneity will give rise to different field-dependences 
of the magnetoresistance. 

We will model the random magnetic field by the following expression: 

B(r)=J2b(r- ri ), (1) 

i 

where r^ is the position of the ith impurity (flux tube, lead grain,...) and b the field produced 
by the impurity (including any constant external field). We will assume that the positions 
are randomly distributed over the sample and denote an average over the positions r t by an 



overline. Defining B = B{r) we can write 

B(r) = B + 5B(r), (2) 

where 5B has zero average. 

III. THE BOLTZMANN EQUATION 

The standard derivation of the Boltzmann equation from the full quantum mechanical 
transport equation relies on a separation of length scales:Q The external fields vary slowly in 
comparison with 1/kp, while the rapidly varying impurity potentials are included in a self- 
energy, giving rise to the collision term in the Boltzmann equation. A previous attempts 
to include random magnetic fields in the Boltzmann equation has viewed the magnetic 
'impurities' as (asymmetric) scatterers, i. e. has included them in the scattering term on the 



right hand side. This, of course, can always be done, even when the correlation length a of 
the modulation (the London length or the size of the lead grains) is much greater than 1/kp. 
However, when (as is the case in the experiments) a ^> 1//c_f, the modulation SB can just 
as well be treated as an ordinary external field in the driving force term of the Boltzmann 
equation. (By contrast it would be inconsistent to put 5B on the left hand side when 
a ~ 1//cf-) In the present paper we will take the latter point of view, i. e. put the random 
magnetic field on the left hand side of the Boltzmann equation. For the ordinary impurity 
scattering we will assume the relaxation time approximation with a constant scattering time 



The Boltzmann equation for the distribution function f(r, v) is then 



Or 



where the Lorentz force ma(r) is 
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a(r) = (E + vxB(r)), 
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We calculate to linear order in the external electric field E, i. e. replace —{e/m)E--^- — ► —eE- 
v dWi wnere f° 1S the Fermi-Dirac distribution. Introducing polar co-ordinates v,(f> for the 
velocity, at T = (the only case considered here) v only enters through -^ = —5{E — Ep) = 
— (mv F )~ 1 5(v — Vp) and can be put equal to Vp. Writing /(r, v, <f>) — f° + g(r, (j>)S(v — Vp) 
we get the following equation for g: 
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Here Uq = eB^jm and 5u(r) = eSB/m. With the following definitions (the i's are for 



convenience : 
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we can write Eq. 



v siiKp ; 
as an operator equation 

Dg = (D + W)g = X - 



The strategy is now to find the Green's function G for D. Then we get 
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where n is the electron density. 
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IV. EIGENFUNCTIONS AND GREEN'S FUNCTION FOR D 

It is straightforward to see that D has the following complete set of eigenf unctions and 
corresponding eigenvalues: 

1 
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Here A is the area of the sample and we have defined 

R{6) = r c 



sine) 



(15) 
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where r c = vf/ujq is the average cyclotron radius. 
The Green's function is then given by 

G (r A y^) = j: ^ r ^ff^ (16) 

- eW-*V<*>r 8 ( r _ r > _ R (0) + ^(0')) ( 17 ) 



cj e 27r /^ T - 1 

(see appendix [A] for details). We use the notation [9] to denote 6 (mod 2ir). Note that G 
is only a function of r — r' because of the translational invariance of D . Inserting Eq. ( |T?D 
in Eqs. (|TT1) and (|l^) for the current, the ordinary Drude formula is recovered. 

V. CALCULATION OF THE RESISTIVITY 

We first have to consider how to perform the average over the random magnetic field. 
Since the distribution function g is the physically relevant quantity, it is g we have to average. 
From Eq. ([11]) we see that this means we have to find the averaged Green's function (%, 
of course, only depends on the external electric field) for D. To do this we consider the 
expansion of D~ l : 

D- 1 = (D + W)- 1 = D G l - D G l WD Q l + D G l WD a l WD Q l . (18) 

We now average Eq. (|H|) term by term. This is formally similar to the quantum mechanical 
treatment of ordinary impurity scattering. Like in quantum mechanics we can organize the 
terms into diagrams as shown in figure |l]. Summing the appropriate geometrical series we 
can rewrite the equation for the averaged Green's function in terms of a self-energy: 



D-i = (D + E)-\ (19) 

where £ is the sum of all irreducible diagrams, i. e. diagrams that do not fall apart when 
an internal line is cut, as shown in figure [| 

We shall only go to second order in W, i. e. only keep the first diagram in figure 0. The 
higher order diagrams (including the maximally crossed ones, of which the first is shown in 

6 



the figure) can in principle be included; they are all finite and we do not expect them to 
change our conclusions (as we have explicitly verified for the second diagram in figure 0).£l 
There are no divergences corresponding to 'weak localization' in this purely semiclassical 
calculation. 

We thus truncate the self-energy and write 



D- 1 = (D - WDo l W)-\ (20) 



The 'self-energy' WD l W is calculated as follows. Using that W is Hermitian we get: 



(nk\WD^W\rik') = / drd(j)dr'd(i)'(Wip nk (r } (i)))*G (r,(i);r' 1 (i)')W^ nlkl (r' , 0') 



dr d<fi dr' d(f>' 8uj(r)5uj(r') 
x[n + k- r(0)] K + k' ■ r(0')]<Mr, 0)*Vvfc'(r', <f>')G (r, 0; r', 0')- (21) 



Here r(0) = r c (cos0,sin0). Now the averaging gives (see appendix |B|) 5u(r)5uj(r') = 
f(r — r'), where the correlation function / (which actually only depends on \r — r'\) depends 
on the nature of the random magnetic field modulation. The averaging restores translational 
invariance, meaning that the self-energy (^l|) (and therefore D~ x ) is diagonal in k. This is 
also seen explicitly by changing integration variables to r\ = r — r' and r2 = r + r' . The 
only ^-dependence of the integrand is then through the factor exp(— i(k — k') ■ T2/2), giving 
a dk,k' upon integration. Therefore G only depends on r — r'. From Eq. (pTTj) we see that 
only (A times) the k = component is needed for the current. We can therefore put k = 
in the self-energy (|2"T|). Inserting the expression ( |TTD for G we can use the delta function to 
do the r*i-integration, giving 
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The integrand is periodic in 0', so we can shift the limits on the ^'-integration to and 
+ 27r. Changing variables to 9 = 0' — and using 

|i?(0)- J R(0')|=2r c |sin^^|, (23) 



the ^-integration gives a 8 n , n i (which was obvious from the start: when k = the matrix 
element is also rotationally invariant.) The final result for the self-energy is 

i 1 r 27T 

E n = (raOlWZ^WlnO) = =-, n 2 f(2r c sm-)e e/uJ ° T e- me d9. (24) 

Now the left hand side of Eq. fl2"0|) is diagonal and we get 

(nO\W\n'0) = 5n f * • (25) 

For the Green's function we get 

/ drG(r, 0; r', <f>') = — £ —L— — . (26) 



It is straightforward to do the remaining angular integrals in Eq. (|lTD and (|T|) to get the 
current density. The result for the resistivity tensor is 
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where we have used S_ n = — S* and defined the renormalised quantities 

uj = ujq - Re Ei (28) 

4 = --ImEi. (29) 

r r 

We see that the change in p is directly related to E x : 

^^ = -rImS 1 (30) 

^^ = _ w - 1 ReE 1 . (31) 



VI. APPLICATIONS 

In ref. |5] we have applied the above theory to the case where the field modulation is 
caused by the Meissner flux expulsion from deposited grains of Pb. Furthermore we also 
showed that the theory could be applied to the experiments by Geim et al.a where the field 



modulation was caused by flux tubes. This establishes that the increase in the resistivity 
even in this case has an essentially classical origin. 

The general analysis of the magnetoresistance can be quite complex owing to the three 
characteristic length scales in the problem: the mean free path l T = vft, the magnetic field 
correlation length /, and the average cyclotron radius r c = vp/uj (when the magnetic field 
is such that quantum effects become important, a fourth length scale, 1/&f, enters). In 
this section we shall confine ourselves to considering the following correlation function for 
different parameter values: 

f\r)=n l 5u^e-^^\ (32) 

which arises from an assumed gaussian field from a single impurity: 

5uo(r) = 5uje- r2/l2 . (33) 

This correlation function is applicable to the experiment in ref. |5], when the modulation 
was dominated by frozen flux in the grains (in this case 5uj would depend on the external 
field and go to zero as B — > B c ), or the case when the modulation is caused by permanent 
magnets. 

We can find an explicit result for the change in the longitudinal resistance for Bq = by 
expanding the integral in the self-energy (|24l) (only 9 ~ 27r contributes). We get 



Afe - (B ° = 0) = (&r)V? fe-(« V -V S . (34) 

PxxO 2 JO 

The overall shape of the magnetoresistance is determined by the ratio x = l T /l. When 
x < 2 there is a maximum in p xx at B = 0, i. e. a negative magnetoresistance. As 
x is increased B = is still a local maximum, but the magnetoresistance develops an 
intermediate minimum and has a maximum for B > 0. The location of the maximum is 
approximately given by the condition that the average cyclotron radius should equal I. When 
x ^ 4, Bo = becomes a local minimum. In fig. |3| we show a plot of the relative change 
of the longitudinal resistance for x varying between 1 and 100. Fig. shows the relative 



change of the Hall resistance for the same range of x'es. By varying the electron density, 
the mobility, and the size I of the impurities, the different types of behaviour might be seen 
experimentally. We note furthermore that it is possible to have a very large enhancement of 
the zero-field resistance in a random magnetic field (for reasonable parameters we can easily 
get a 20-fold increase). 

The problem of the (quantum mechanical) localization of a particle in a random magnetic 
field has received a lot of attention.^ As we have seen we can get a negative magnetoresis- 
tance around Bq = 0, reminiscent of weak localization. This enhancement of the resistance 
could be interpreted as a classical sign of 'localization' in the random magnetic field. How- 
ever, varying the parameters it is also possible to get a local minimum around B = 0, so 
there is probably no simple connection to the quantum mechanical theory of localization 
(although we always get an enhancement of the resistance). 

We end this section by some speculations regarding the possibility of applying our theory 
to the Quantum Hall v = 1/2-state. As shown in the seminal paper by Halperin et alJiil it 
is possible to make a Chern-Simons gauge transformation and transform the v = 1/2-state 
to a system of fermions (socalled composite fermions, which can be viewed as an electron 
with two flux-tubes attached to it) in zero average magnetic field. The composite fermions 
interact through a fictitious Chern-Simons magnetic field given by 

5b cs =2%5p, (35) 

where $o — h/e is the flux quantum and Sp is the deviation of the electron density from the 
average density n. In a static approximation 5p is nonzero simply because of the electrostatic 
coupling to the impurities. We then have a system of fermions moving in a static, random 
magnetic field, i. e. the system considered above. If the impurities are fully ionized we would 
have rii = n, i. e. a much higher impurity density than in the grain experiments, and I is 
typically much smaller than before (of the order of 500 A). In our model this will give rise 
to a broad minimum in p xx around B = (that is around -B1/2) for our simple model 
correlation function (see fig. |^). 
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In fact, we are here dealing with the limit x ^> 1, in which case the self-energy ( p4j ) can 
be evaluated analytically for moderate fields. We get 



3 



^=^Vn/^coth(^ (36) 

This is consistent with the behaviour seen experimentallyJi3TL3 We also note that typical 
values of p u =i/2/po are in the range 25-500. To do quantitative comparisons with experiment 
there is an additional complication: Since the impurities give rise both to the magnetic field 
and the impurity scattering it is inconsistent to put them on different sides of the Boltzmann 
equation, when in reality they are strongly correlated. Instead the potential scattering should 
be brought to the LHS of the Boltzmann equation as a random electric field. This does not 



contradict the discussion in the beginning of section |T|, since kpd ~ 15 ^> 1 in a typical 
QHE sample.tll (d ~ the distance from the 2DEG to the impurity layer, sets the distance 
over which the impurity potential varies appreciably). Calculations are currently underway, 
addressing this issue. 

VII. CONCLUSION 

We have shown how to calculate the magnetoresistance of a 2DEG in a random magnetic 
field in the semiclassical approximation. The magnetic field was included in the driving 
force term in the Boltzmann equation. Furthermore we have applied our results to different 
experiments and shown how the different types of random magnetic fields gives rise to the 
differences in the magnetoresistance. Finally we speculated on the possible relevance of our 
theory to the v = 1/2 Quantum Hall state. 

We would like to thank Mads Nielsen and Kasper Juel Eriksen for discussions. 



APPENDIX A: GREEN'S FUNCTION FOR D, 
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Inserting the expression (|I3"D for the eigenfunctions into Eq. flTB]) we get 
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G (rA;r>,4>') = ^E rv ^E^ (r "^ W)) - ( A1 ) 

The fc-sum just gives a delta function. The sum over n involves calculating 



27r „ n + ia 

1 °° r 00 1 

= — Y e i(2 ™ + ^^^rfx, (A2) 

where we have used the Poisson summation formula. The integrand has a simple pole in 
x = —ia. If 27rm + > 0, we close the contour in the upper half plane and get 0. If 
2nm + (p < 0, we close in the lower and pick up a contribution — 2iiie^ 2nm+ ^ a from the pole: 

h((f>) = -i J2 e (27rm+ * )a . (A3) 

2ivra< — <j> 

Since h is periodic in <fi we can choose the argument to lie between and 2n and get 



-i 



h{4>) = -ie a[ ^ E e 
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= -^ m Y—^- (A4) 

With a = 1/uqT we then get Eq. (0). 

APPENDIX B: CALCULATION OF THE CORRELATION FUNCTION 



We have to calculate 5uj(r)5uj(r') = (e 2 /m 2 )5B(r)5B(r'). The averaging consists of 
integrating over the iV positions (assumed independent) of the impurities. We get 



5B(r)5B(r') = B(r)B(r f ) - B 2 



= Y /b (r-r l )b(r'-r J )-B 2 

1 r N 

= -JW E / II ** Kr - r i )6(r' - r 3 ) - B 2 . (Bl) 

71 i,j fc=l 

The N(N - 1) terms with i =£ j are equal to A N ~ 2 (J dRb(R)) 2 = A N B 2 /N 2 (since B = 
B = No), and the N terms with i = j are equal to A N ~ l J dRb{r — R)b(r' — R). Dropping 
the term of order 1/N we get 
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8B{r)8B{r') = m f dRb(r - R)b(r' - R), (B2) 

where rij = N/A is the density of impurities. The correlation function thus becomes 

f(r)=m f dRb(r + R)b(R). (B3) 

It is only a function of r. Note that it does not matter if we write Sb for b in (|B3|) , as b and 
5b only differ by a constant of order 1/N. 
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FIGURES 
FIG. 1. The perturbation series for the full propagator (Green's function) G. The unperturbed 

propagator Go is denoted by a thin arrow. A cross denotes W, and the dashed lines an impurity 

average. 

FIG. 2. The (irreducible) self-energy X. The calculations in the text approximates T, by the 
first diagram in the series. 

FIG. 3. Relative change in longitudinal resistance of a 2DEG in a random magnetic field with 
the correlation function p^). Parameters as follows: n = 4 10 15 m~ 2 , ni = 10 9 m -2 , db = 0.03 T. 
The curves show (from the top) x = 1, 2, 3, 5, 10, 20, 100. We are varying I and l T such that ll T = Z 2 ,, 
where Iq = 8.1 /im is constant, to keep the curves on the same scale. 

FIG. 4. Relative change in Hall resistance for x varying between 1 and 100, with parameters 
as in fig. [3|. x = 1 corresponds to the lowest curve. 

FIG. 5. The relative change in the longitudinal resistance with parameters as follows: 
n = 8 10 15 m~ 2 , m = 8 10 15 m -2 , I = 500 A, fi = 100 m 2 /Vs, db = 0.03 T. The change in 
the Hall resistance (not shown) is very small, of the order of 10 -3 . 
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